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•RELATIVELY  OPTIMAL  REINFORCEMENT  PATTERNS 
FOR  FIBER-REINFORCED  COMPOSITE  MEMBRANES 

A.  A.  G.  Cooper*  and  E.  M.  Wu** 

Materials  Research  Laboratory 
Department  of  Mechanical  and  Aerospace  Engineering 
Washington  University 
St.  Louis,  Missouri  63130 

I .  Summary 

A  method  is  proposed  to  obtain  relatively  optimal  reinforcements  in  a  fiber- 
reinforced  composite  membrane.  The  method  is  based  on  the  fact  that: 

a.  A  reinforcing  fiber  is  most  efficiently  utilized  when  it  coincides 
with  the  direction  of  maximum  required  stiffness. 

b.  The  directions  of  principal  trajectories  are  not  dependent  on  any  fiber 
reinforcement  as  long  as  the  fibers  coincide  with  those  directions. 

An  optimality  condition  for  the  fiber  reinforcement  is  derived.  The  derivation 
is  an  adaptation  of  one  given  in  the  literature  for  isotropic  homogeneous  materials.  The 
optimality  condition  derived  states  that  the  specific  strain  energy  divided  by  the  speci¬ 
fic  number  of  fibers  is  a  constant,  for  maximum  stiffness  at  a  given  weight,  and  mini¬ 
mum  weight  for  given  stiffness.  Since  optimization  of  the  reinforcement  only  is  con¬ 
sidered  and  optimization  of  the  matrix  is  not  included,  the  resulting  optimum  is  relative. 

Mathematically,  the  optimum  is  weak.  As  an  example,  some  numerical  values 
for  a  boron-epoxy  wedge  in  tension  were  computed. 
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3.  Notation. 

(Notation  used  in  appendix,  not  included) 
(7)  equation  (7) 

[7]  Reference  number  171 
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Pj  external  load;  f  =  1,2, _ 

rj  radius;  figure  1;  1  =  1,2. 

5  compliance  coefficient;  I,  i  =  l,2f6. 

•I 

U  specific  strain  energy  (1 8.2.) 

U.  component  of  U,  dependent  on  a.  (f  8.2.) 

V  volume 

WD  work  done  by  external  loads 

D 

a  wedge  angle;  figure  1. 

6  variation 

8j  displacement  of  loadpofnt  d F  Pj 

strain;  I  =  1,2,6. 

oj  stress;  I  =  1,2,6. 

4.  Nomenclature; 

FLAG:  Fiber  LAylng  Gadget. 

FRC:  Fiber  Reinforced  Composite 
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RELATIVELY  OPTIMAL  REINFORCEMENT  PATTERNS 
FOR  FIBER-REINFORCED  COMPOSITE  MEMBRANES 

A.  A.  G.  Cooper 
Materials  Research  Laboratory 
Washington  University 
St.  Louis,  Missouri  63130 


5.  Introduction 


Fiber-reinforced  composites  (FRC's)  can  be  of  great  advantage  in  structures  or 
structural  parts  requiring  a  material  with  distinct  directional  characteristics.  Pressure 
vessels,  for  instance,  are  such  structures,  and  fiber-wound  pressure  vessels  are  examples 
of  successful  utilization  of  FRC's. 


The  directional  load  or  stiffness  requiremenrs  are  not  always  constant  throughout 
the  structure;  quite  often  they  change  rapidly,  i.e.  the  directions  and  values  vary  con- 

siderably  over  a  short  distance.  A  well-known  example  is  a  strip  with  a  hole  under 
tension. 


FRC's  are  most  commonly  applied  in  the  form  of  cloth  or  tape  made  of  continu¬ 
ous  straight  fibers.  Therefore,  in  order  to  follow  changes  in  directions,  a  laminate  has  to 
be  made  consisting  of  several  layers  with  different  orientations.  It  is  obvious  that  the 
material  is  not  used  very  efficiently  in  these  cases  of  rapidly  changing  directions,  the 
more  so  since  generally  the  directional  loads  and/or  stiffness  requirements  vary  simul¬ 
taneously  and  significantly  as  well.  Hence  it  appears  to  be  worthwhile  to  look  into  the 
possibility  of  a  more  satisfactory  method  to  deal  with  this  kind  of  problems. 

Simply  stated,  the  problem  is  how  to  use  as  little  material  as  possible  to  provide 
as  high  a  stiffness  as  possible  in  a  given  direction,  and  clearly  this  is  an  optimization 
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problem. 

In  this  report  an  account  is  given  of  initial  work  d>ne  on  a  method  to  optimize 
the  reinforcement.  It  is  based  on  two  observations: 

a.  For  a  uni  directionally  reinforced  FRC,  the  direction  of  maximum  stiff¬ 
ness  coincides  with  the  fiber  direction. 

b.  The  principal  stress  (and  strain)  directions  in  an  isotropic  body  do  not 
change  when  the  isotropic  body  is  made  specially  orthotropic;  i.e. 
when  the  axes  of  orthotropy  coincide  with  the  principal  stress  (strain) 
directions. 

By  laying  fibers  along  principal  trajectories  determined  for  the  isotropic  case, 
reinforcements  are  obtained  which  are  (relatively)  optimal  under  certain  conditions  for 
the  stress  anchor  strain  fields.  This  is  the  essence  of  the  method  discussed  in  Section 
of  this  report. 

Optimum  designs  are  often  impracticable  because  they  are  too  complicated  and 
hence  too  expensive  to  fabricate.  This  is  not  likely  to  happen  to  FRC  designs  optimized 
as  outlined  in  Section  8,  for  the  following  reasons: 

a.  The  price  of  a  reinforcing  material  like  boron  or  graphite  is  so  high  [35]  that  it 
contributes  substantially  to  the  cost  of  the  structure.  Any  weight  reduction  causes 
a  similar  cost  reduction  of  the  reinforcing  matejjal  and  hence  also  an  overall  cost 
reduction  for  the  structure. 

b.  One  of  the  main  factors  contributing  to  material  cost  is  the  precision  required  for 
fiber  alignment.  High  precision  is  also  required  when  the  material  is  applied  to 
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the  structure.  For  the  method  discussed  in  Section  8,  the  material  is  used  in  the 
form  of  a  single  continuous  wire.  Hence  the  need  for  precision  during  fabrication 
has  been  eliminated,  and  of  the  two  precision  requirements,  only  the  one  for  appli¬ 
cation  remains.  Not  only  will  this  result  in  lower  cost  but  it  will  be  beneficial  also 
for  the  overall  precision  finally  achieved. 

c.  The  percentage  of  material  waste  will  be  very  low  for  production  versions  of  the 
fiber-laying  device  described  in  Section  9,  since  the  length  of  the  fiber  needed  can 
be  computed  accurately.  (Production  versions  of  the  FLAG  will  be  equipped  with 
a  fiber  cutting  device.  Therefore,  no  fiber  will  be  laid  along  return  loops.) 

In  order  to  put  the  method  in  a  proper  perspective,  some  optimization  theories 
and  techniques  used  in  the  field  of  FRC's  as  well  as  a  few  other  fields  are  briefly 
surveyed. 

Finally,  a  description  is  given  of  an  experimental  Fiber-LAying  Gadget  (FLAG) 
which  is  being  developed  to  lay  s  single  continuous  fiber  along  a  curved  trajectory  in  a 
flat  plane.  This  device  will  eventually  be  capable  of  laying  fibers  in  predetermined 
patterns,  and  then  test  specimens  can  be  made  to  check  theoretical  results  with  experi¬ 
mental  data. 

6.  Structural  Application  Methods  for  FRC's 

In  this  section  a  brief  survey  is  given  of  methods  used  to  obtain  optimal  struc¬ 
tural  applications  of  composites  reinforced  with  continuous  fibers.  Of  these  methods, 
drawbacks  are  indicated  which  have  been  aleviated  by  the  method  discussed  in  Section  8. 
Extensive  bibliographies  covering  each  of  these  methods  can  be  found  elsewhere,  see  for 
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Instance  Reference  1. 

6.1 «  Lamination 

Continuous  fibers  are  often  worked  into  tapes  or  woven  fabrics.  In  a  laminate 
made  of  layers  of  tapes  anchor  fabrics,  the  number  and  the  orientation  of  the  layers  are 
variables  which  allow  the  structural  designer  to  more  or  less  match  the  material  proper¬ 
ties  and  the  structural  requirements. 

Methods  hove  been  developed  to  determine  rationally  the  number  end  orienta- 
tion  of  the  layer.  ([2,  3,  4]  end  [5],  Dlv.  6.  Propulsion  Systems).  Apert  from  the  feet 
thot  sometimes  these  method,  reportedly  yield  questioneble  results,  they  determine  the 
varlebles  in  one  single  point  only.  In  appllcetlon.  where  the  stress  end/or  strain  gradient 
and  direction  do  not  vary  considerably  throughout  the  structure,  them  method,  will  pro¬ 
vide  useful  infomretlon.  However,  in  structure,  with  stress  concentration,  they  become 
imprecticel  since  they  require  epplic.tlon  of  the  material  in  quantities  and  form,  that 
are  impracticable.  The  designer  therefore  is  forced  to  compensete  by  using  more  mat  >- 
rial  and  thus  decreasing  its  efficiency. 

As  a  result,  the  superior  directional  properties  are  not  used  to  full  advantage  in 
many  laminated  structures. 

6.2.  Filament  Winding 

Fiber-wound  structures  account  for  some  of  the  most  successful  applications  of 
FRC's.  This  is  mainly  due  to  the  uniformity  of  the  stress-and-strain  fields  in  structures 
particularly  suited  for  fiber  winding,  e.g.  pressure  vessels. 
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Fiber-winding  techniques,  to  a  certain  extent,  allow  for  aligning  fibers  with 
principal  directions  and  controlling  the  number  of  fibers  in  a  given  direction.  In  some 
cases  geometrical  disturbances  like  openings  in  vessels  can  be  incorporated  in  the  design 
in  such  a  manner  that  the  design  is  optimal  [1] . 

Stress  analysis  of  fiber-wound  structures  is  based  either  on  netting  analysis  or  on 
a  continuum-mechanics  approach  [6].  In  netting  analysis,  most  frequently  used  in  the 
early  days  of  filament  winding,  the  fibers  are  assumed  to  be  the  load-carrying  consti¬ 
tuents  of  the  composite,  with  the  matrix  not  playing  a  load-carrying  role  at  all.  Netting 
analysis  simplifies  the  stress  analysis  of  fiber-wound  structures  considerably  and  in  many 
papers  its  application  has  been  investigated,  e.g.  [7-11] .  Studies  of  some  special  struc¬ 
tures,  not  being  fiber-wound  structures,  have  also  been  based  on  netting  analysis  [12]. 

When  based  on  continuum  mechanics,  the  analysis  of  fiber-wound  structures 
becomes  complicated  [6].  Since  the  structures  under  consideration  have  one  or  more 
axes  of  rotatory  symmetry  and  also  because  the  stress  and  strain  fields  are  un  form,  the 
complicated  analyses  are  still  tractable. 

6.3.  Reinforced  Concrete 

Reinforced  concrete  can  be  classified  as  a  fiber-reinforced  composite.  The 
functions  of  reinforcement  and  matrix,  as  assumed  in  analyses  of  reinforced  concrete 
structures,  differ  from  those  in  "real"  FRC's;  the  reinforcement  is  assumed  to  carry 
tensile  loads  while  the  matrix  carries  compressive  loads. 

The  reason  for  paying  attention  to  the  field  of  reinforced  concrete  is  the  fact 
that  in  that  field,  too,  there  is  a  need  for  optimization  of  the  reinforcement.  And  some 
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interesting  methods  have  been  developed,  both  theoretical  [13]  and  experimental  [14]. 

In  [14]  a  discussion  has  been  given  on  the  determination  of  the  reinforcement  pattern 

in  a  reinforced  concrete  plate  of  irregular  form  by  employing  Ligtenberg's  reflection- 

moire  method.  With  such  an  experimental  method,  the  method  discussed  in  Section  8 

« 

could  still  be  used  in  cases  where  analytical  and  numerical  analyses  are  impractical. 

7.  Optimum  Design 

There  is  a  distinct  difference  between  "structural  analysis"  and  "structural 
design."  Analysis  is  the  study  of  the  behavior,  i.e.  stresses,  deflections,  etc.  of  a  struc¬ 
ture  of  given  form,  dimensions,  and  material  constants.  Design,  or  direct  design  as  it 
sometimes  is  called,  is  concerned  with  structures  in  which  certain  parameters,  e.g. 
plate  thickness,  are  to  be  determined.  These  parameters  are  determined  such  that  the 
structure  not  only  satisfies  conditions  of  equilibrium  and  compatibility  but  also  meets 
one  or  more  additional  requirements.  For  example,  the  requirement  that  the  stiffness 
be  maximum  for  a  given  weight  or  the  weight  be  minimum  for  a  given  stiffness.  The 
structure  resulting  from  such  a  structural  design  procedure  is  an  optimum  (or  optimal) 
design. 

On  the  subject  of  optimal  structural  design,  a  most  extensive  literature  now 
exists  and  quite  a  few  survey  papers  cover  parts  of  the  subject.  Barnett's  paper  [15] 
is  a  good  introduction  (and  a  very  readable  one  at  that)  to  some  optimum  design  topics 
related  to  the  method  discussed  in  Section  8.  Also,  a  survey  article  [16]  is  mainly 
devoted  to  reinforce  J  concrete  optimisation. 

The  following  paragraphs  in  this  section  refer  to  fields  in  optimum  design  theory 
which  are  related  to  the  problem  of  optimization  of  FRC's. 


7.1.  Plastic  Design 


Plastic  design  deals  with  structures  In  which  plastic  deformation  Is  allowed  to 
occur.  Much  progress  has  been  made  in  developing  the  theory  and  many  important 
results  have  been  established  concerning  optimality  conditions. 

In  [17],  minimum  weight  plastic  design  has  been  related  to  uniform  strength 
elastic  design,  based  on  the  assumption  that  the  yield  function  is  directly  proportional 
to  the  specific  strain  energy.  One  of  the  important  theorems  in  the  theory  of  plastic 
design  states  that  a  design  is  a  minimum  weight  design  if  the  rate  of  dissipation  per 
unit  volume  is  constant.  Analog  to  this,  for  an  elastic  design  to  be  an  optimum  design 
the  specific  strain  energy  is  constant  [17]. 

7.2.  Elastic  Design 

Derivation.  of  optimality  conditlom  for  alaitlc  .ttuchir..  at.  ofrtn  bawd  on  th. 
principle  of  minimum  pot.nti.l  energy.  Thi.  prlnctpl.  .t.tn  (m.  for  Imtanc.  [18], 
pag.  171)  that  th.  dllplac.rn.nt  which  Mti.fi..  th.  differential  equation,  of  .ppilibrium, 
a.  well  a.  th.  condition,  at  th.  bounding  surface,  yield,  a  smaller  value  for  th.  potential 
energy  of  deformation  than  any  other  dl.pl.cm.nt,  which  .ati.fi..  the  Mm.  condition. 

at  the  bounding  surface. 

Usually,  only  one  type  of  structure  is  treated  at  a  time,  e.g.  sandwich  beam, 
column,  etc.,  and  a  linear  relationship  between  stiffness  and  weight  is  assumed  [19] . 
Furthermore,  in  most  cases  only  one  design  parameter,  like  the  thickness  of  a  plate,  is 

independently  variable. 
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A  generally  applicable  optimality  condition  for  maximum  stiffness  designs  has 
been  derived  by  Huang  [20] .  His  very  simple  derivation,  also  based  on  the  principle  of 
minimum  potential  energy,  leads  to  a  condition  for  optimality.  This  derivation  is  used 
in  Section  8  to  establish  an  optimality  condition  for  the  reinforcement  in  FRC  membranes. 

Since  an  optimization  problem  can  be  formulated  as  an  isoperimetric  problem 
[21,  22],  variational  calculus  is  frequently  used  in  derivations  of  optimality  conditions. 
Once  an  optimality  conH? non  has  been  established,  for  instance  to  obtain  maximum  stiff¬ 
ness  for  a  given  weight,  another  optimality  problem  with  the  same  condition  can  easily 
be  formulated  to  obtain  minimum  weight  for  a  given  stiffness.  This  is  due  to  what  is 
called  t..u  reciprocity  of  isoperimetric  problems  [21];  in  the  theory  of  mathematical 
programming  it  is  known  as  duality  (see  for  instance  [23]).  This  approach  has  beon  fol¬ 
lowed  in  [24] . 

7.3.  Maxwell -Michel  I  Structures 

A  theory  of  optimum  frame  structures  based  on  work  by  Maxwell  [25],  has  been 
developed  by  Michel  I  [26].  According  to  this  theory,  the  bars  of  an  optimum  frame 
coincide  with  the  principal  strain  trajectories  of  a  displacement  field  involving  constant 
principal  strains  the  value  of  which  is  determined  by  means  of  Hooke's  Law  by  the 
allowable  stress  and  the  E-modulus  of  the  bar  material.  The  principal  strain  lines  can 
be  shown  to  be  given  by  the  same  equations  as  the  slip  lines  in  plane  flow  of  a  perfectly 
plastic  solid  [27,  28] .  For  the  story  of  the  revival  of  the  interest  in  Michel  I  structures, 
see  [15,  27] . 

In  [29]  it  has  been  suggested  that  this  theory  can  be  applied  to  FRC's.  A  lami¬ 
nate  of  varying  number  of  layers  has  to  be  used  to  approximate  the  change  in  fiber 
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diameter  required  by  the  theory.  No  justification  has  been  given  for  this  application 
of  the  theory. 

Michell  -*uctures  are  statically  determinate  and  can  therefore  be  dimensioned 
such  that  the  designs  rre  uniform  strength  designs  md  optimum  designs  as  well.  Gener¬ 
ally,  FRC  structures,  e.g.  plates  under  a  biaxial  state  of  stress,  are  statically  indetermin¬ 
ate  and  consequently  applicability  of  the  theory  to  FRC  structures  is  far  from  being 
self-evident. 

7.4.  Systematic  Structural  Synthesis 

The  methods  of  optimization  mentioned  in  the  three  preceding  paragraphs  are 
analytical.  However,  much  work  of  a  numerical  nature  has  also  been  done.  By  means 
of  numerical  methods,  cases  can  be  investigated  which  are  too  complicated  to  be  treated 
analytically,  e.g.  highly  indeterminate  structures. 

Systematic  structural  synthesis  is  the  procedure  leading  to  an  optimum  design  by 
iteration  of  an  analysis-redesign  cycle.  For  the  iteration,  quite  a  few  numerical  methods 
are  known  [1] .  Ths  primary  requirement  to  be  met  bv  such  methods  is  a  rapid  conver¬ 
gence  to  an  absolute  optimum. 

Structural  synthesis  is  basically  the  same  trial  and  error  design  procedure  as  has 
always  been  used  by  structural  designers.  It  differs  in  that  the  procedure  has  been  ration¬ 
alized  in  the  sense  that  in  the  iteration  process  every  next  step  yields  a  design  that  is 
better  than  or  at  least  as  good  as  the  previous  one.  This  is  achieved  by  employing  mathe¬ 
matical  methods  developed  in  the  theory  of  mathematical  programming  [23,  301.  Because 
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„  *.  enormous  isoi,  of  numerical  operation,  involved, 

,f  impute,. _ *  *™n  then,  particularly  whan  noo-Un-r  pro».mmlnfl 

method,  are  applied,  only  a  limited  n-ber  of  design  variable,  can  be  handled  to  pre¬ 
vent  ,h.  problem  from  biomlng  too  implicated  to  formulate  [311. 

[321,  systematic  synthesis  ha.  bin  discoid  and  lll^.t.d  by  an  example 
problem.  An  Important  reH.lt  reported  In  the  paper  I.  the  obirvation  that  the  minimum 
optimum  detlgn  fo,  a  ...tlcally  Indetiminat.  .tructur.  I.  no,  nieHarl.y  a  fully 

stressed  design. 

The  significance  of  statical  determinccy  or  statical  Indeterminacy  and  fully 

i  v  ‘eatinn  was  not  very  clear  in  Hie  past  [15],  but  has 
stressed  design  in  structural  optimization  was  n  y 

since  bin  ini.tlg.ted  thorourfsly  [16,  3,1=  A  statically  lnde«rmin.«  structure  1.  • 
relatively  optima,  design,  which  „  no,  fully  .Hiid  and  no,  stiff.,  for  die  sum.  weight 
nor  lighter  for  the  same  stlffnit  than  a  corrwponding  st.Hc.lly  deWminate  opt,  urn 

design. 

Sometimes  a  statically  indatwmlnat.  structure  can  become  fully  .,«»d  if  it 
h„  bin  prestressed  [331.  This  fully  strewed  d«ign  is  llfr  *■"  *“ 
obtained  without  prestressing. 

|f  ,h6  domain  of  fislbl.  Kslution.  I.  no,  convex  or  the  mathinatical  program¬ 
ming  technique,  do  n6t.  In  gieral,  yield  abilu^ly  optimal  solutions,  then  other  method, 

have  to  be  used  [34] . 


From  ,h.  literature  survey,  some  conclusion,  can  be  drawn  which  indicate  pos¬ 
sible  direction,  of  investigation  of  the  optimizaHon  of  reinforcement  in  FRCS. 
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St.tlc.lly  indeterminate  optimal  design,  .re.  In  general,  r.l.ti»ely  optim.l  in 
the  HM  that  an  equivalent  but  rt.tlc.lly  determinate  d«ign  can  be  found  which  i. 
abrt.lut.ly  optimal.  An  orthogonally  stiffened  plat,  in  some  Instance,  can  be  considered 
a,  a  statically  dahmoinate  structure  since  tb.  mo  equilibrium  condition,  cm  be  taken 
care  of  by  tbe  mo  set.  of  appropriately  dimensioned  stiff—  [361.  Unfortunately,  this 
appears  to  be  impossible  for  an  FRC  membrane  because  elastic  inaction  bemeen  two 
mutually  orthogonal  unidirectional  fiber  system,  in  FRC,  would  couple  the  solution,  for 
each  of  the  mo  equilibrium  conditions.  For  the  same  re  »on,  Michell's  theory  doe.  not 
^  generally  applicable  to  FRC  , rtuctures  reinforced  with  fiber,  of  constant  diameter. 
Them  observations  indie...  that  probably  in  only  a  few  special  case,  reinforcements  pat¬ 
terns  can  be  obtained  which  are  absolutely  optimal. 

The  use  of  numerical  methods,  particularly  non-linear  programming  moth  k, 
complicates  matters  considerably.  Therefore,  it  appears  to  be  expedient,  at  least  ini¬ 
tially,  to  investigate  the  optimisation  proelem  analytically  rather  then  numerically. 

An  unsuccessful  attempt  was  mad.  to  find  way.  to  apply  the  rault.  of  f  171  to 
FRC1,.  The  elastic  strength  concept  a,  mentioned  in  [17]  retire,  the  stnrctur.  to  be 
statically  determinate.  This  restriction  seventy  limits  the  applic  r  ?  <*  >"  op"' 
reality  condition  derived  from  [171.  Therefore,  this  approach  ha,  not  bee,  pursued. 
Moreover,  completely  satisfactory  strength  theories  for  composite,  have  not  yet  be 
established.  Any  optimization  thexy  based  on  a  strength  theory  could  be  expect,  to 
be  plagued  by  such  nuisances  as  uon-inv.riu.cy  with  r«pect  to  coordinate  tr.nrfom.a- 
Hon  or  ambiguity  of  principal  directions. 
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Numerical  methods  are  virtually  Inevitable  when  a  structure  Is  to  be  subfected 
to  many  different  loading  conditions.  Such  a  cate,  however,  will  be  considered  to  be 
beyond  the  scope  of  this  report,  which  is  primarily  meant  as  an  attempt  to  grasp  the 
fundamentals  of  reinforcement  optimization. 

8.  Optimization  of  Fiber  Patterns  in  FRC  Membranes 

In  this  section,  after  some  preliminary  remarks,  an  optimality  condition  is 
derived  for  FRC  membranes.  The  derivation  is  a  direct  adaptation  of  the  one  given  by 
Huang  [20]  for  homogeneous  isotropic  materials.  The  term  membrane  refers  to  struc¬ 
tures  with  in-plane  stress  resultants  only;  it  implies  that  coupling  effects  are  not  to  occur. 
Hence,  only  balanced  laminates  are  considered. 

In  the  third  paragraph  of  this  section  a  few  comments  are  made  on  an  illustra¬ 
tive  example  of  seemingly  trivial  form:  a  wedge  in  tension.  Details  of  the  analysis  of 
this  wedge  ap pear  in  the  Appendix. 

8.1.  Preliminaries 

a.  From  experiments  as  well  as  from  theory  it  is  evident  that  in  unidirectional ly  rein- 
force  FRC's  the  direction  of  maximum  stiffness  coincides  with  the  fiber  direction. 
Hence,  if  the  FRC  is  required  to  provide  a  certain  stiffness  in  a  specific  direction, 
the  fiber  material  is  most  efficiently  utilized  when  the  fibers  coincide  with  that 
direction  [35]. 

b.  Directional  reinforcement  does  not  affect  principal  stress  or  strain  directions  as  long 
as  the  directions  of  the  reinforcement  coincide  with  the  principal  directions.  That 
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is,  the  principal  directions  in  an  isotropic  body  remain  unchanged  when  fibers  are 
laid  along  those  directions.  As  in  the  isotropic  case,  the  principal  stress  trajectories 
continue  to  coincide  with  the  principal  strain  trajectories.  This  follows  from  Hooke's 
law: 

*6  =  S14’l  +  S24"2  +  *66’i 

where  e,  =  shear  strain 
o 

a] '  a 2  =  norma^  s,Tess 

a  =  shear  stress 
6 

S,.,S_,,S,Z  =  compliance  coefficients 
1 6  26  66 

If  the  1-  and  2-directions  are  principal  stress  directions,  then  =  0.  Since  the 
principal  lines  form  an  orthogonal  system,  and  the  fibers  are  laid  along  the  principal 
lines,  the  FRC  is  locally  specially  orthotropic.  Therefore,  S^=  5^  =  0.  Thus, 
e6=  0  regardless  of  the  elastic  constants.  That  is,  regardless  of  the  number  of  fibers 
if  the  fibers  coincide  with  the  principal  stress  lines.  As  =  0  is  a  sufficient  condi¬ 
tion  for  the  strains  to  be  principal  strains,  it  follows  that  principal  stress  lines  and 
principal  strain  lines,  which  coincide  in  the  isotropic  case,  continue  to  coincide  in 
the  case  of  locally  special  orthotropy. 

Observations  a.  and  L  lead  to  the  hypothesis  that  fiber  reinforcements  are  opti¬ 
mal  if  the  fibers  follow  the  trajectories.  And  since  the  trajectories  are  not  changed  by 
the  presence  of  the  fibers,  they  are  the  same  as  if  the  material  were  isotropic.  This  is 
fortunate  since  for  the  determination  of  trajectories  in  isotropic  bodies  analytical, 
numerical  and  experimental  methods  are  known. 
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8.2.  Optimality  Condition  for  FRC  Membranes 

An  Infinitesimally  small  element  of  an  FRC  membrane  is  considered.  The 
reinforcement  is  treated  as  a  specially  orthotropic  medium  in  *he  form  of  two  kinds  of 
layers.  In  one  kind  of  layer  the  fibers  coincide  with  the  1-direction  and  in  the  second 
kind  the  fibers  coincide  with  the  2-direction.  The  1-  and  2-directions  coincide  with 
the  principal  stress  (strain)  directions.  Let  the  two  sets  of  layers  be  represented  in  the 
form  of  two  layers.  (In  an  actual  laminate  the  layers  of  the  two  sets  can  bo  arranged 
such  that  no  coupling  effects  occur.)  The  elastic  constants  for  a  layer  are: 


C  (k)  =  C  a  1 
C11  CE  k 


C  (k)  =  C  a  3 
c12  *“P  k 


'22 


<k>  -  C.  a,  2 


'E  “k 


where  C.^  =  stiffness  coefficient  for  kth  layer  (k=  1,  21 
Cr ,  Cp  =  constants 

ak  =  specific  number  of  fibers  in  kth  layer  (specific  =  per  unit  area) 
n.  =  exponent,  to  be  determined  theoretically  or  experimentally 

In  the  following,  by  strain  energy,  the  strain  energy  accumulated  in  the  reinforcement 
is  meant.  The  strain  energy  per  unit  area  of  a  membrane  of  unit  thickness  is: 

U  =  1/2  (oj  e1  +  o2  «2) 


where  o.  and  e.  are  the  principal  stresses  and  strains. 


(with  the  usual  assumption  of  strain  compatibility  between  the  layers) 


r 

•  0) 

(2) 

(D 

(2) 

=  1/2  [ 

,<cn 

+  C11> 

•i  + 

2<C12 

+  cl2>  *1  *: 

=  1/2 

[(bi  *i 

"l  + 
a1  + 

b2  e3 

n3 

a1  + 

n2 

btVl  ) 

+  (b1  e1  e 

"2  + 
'2  + 

b2*3' 

"3  + 

'2  + 

bl  *2  a2  >] 

=  1/2 

[<2Ui> 

+  (2U2)] 

where  Uj  =  component  of  U,  depending  on  a. 


'l  *2 


0) 


Adapting  Huang's  [20]  analysis,  an  optimality  condition  for  FRC's  is  derived  in  the 
following.  Tiie  external  work  Wg  can  be  considered  as  a  measure  of  stiffness.  Wg  can 


be  written  as  [20] : 


-WB  =  / 1/2  (U,  +  U2)dS  -  R  S. 


where  Pj  =  external  loads 

=  displacement  of  loadpoint  of  P. 
dS  =  surface  area  element 


summation  over  repeated  index 
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Two  designs,  B  and  B,  are  compared  which  are  subjected  to  the  same  loading  and  have 
the  same  stiffness,  i.e.  Wg  =  Wg .  Thus, 


1/2 /K  *lal'  +  b2  *3  *1  3  +  bte2"l2) 

+  (blV22  +  b2e3a23  +  ble2s2')]  dS  '  Pi*i 

.i*  r  ni  n3  n'> 

=  1/2  J  [<b,  5,  5,  +  b2  53  S,  +  b,  52  !,  2) 


+  (b,  S,  S22  +  b2S3i23  +  b!  s2  s2  ’ >J  d  S  ’  p,  5| 


By  the  principle  of  minimum  potential  energy  [18],  the  following  inequality  holds: 


'/2/[bi 


nl  ni  1  - 

e,j  +  ...  +  J  d  S  -  P.  5. 


<  l/2/[b,.,i,'  +  ...  +  b,.2i22]dS  -  P.6, 


Substitution  of  (2)  into  (3)  yields, 


n,  n 


n.  n„ 


n^  n. 


,/2/[b,  •l(il'-,)')  +  b2 '3  (i,  3  '  "l  3)  +  b. '2  (i,  2  '  al  2) 


n9  n9  "q  "q  "1  "i  1 

+  b1  e1  (aj  ~  a2'  +  b2  ®3  $2  ~  a2  ^  +  bl  ®2  ^S2  "  a2  d  S  “  ° 


n„  n 


n,  n. 


Assuming  that  a.  and  a.  do  not  differ  very  much,  i.e.  assuming  that  5  a.  is  very 
small  in, 


22 


a.  =  a.  +  8  a. 

i  :  i 


a.n  then  can  be  expanded  as  follows: 


Sn  =  e.'  +  na.  8  a.  +  terms  of  higher  orders  of  8  a.,  which  will  be  neglected 
i  i  i  i  i 


_n  n  n-1  c 

a.  =  a.  +  n  a.  8  a. 

!  I  II 


With  (5),  (4)  becomes: 


n,  -1 


n„  - 1 


n„  - 1 


'/2f  [<bl  *1  "l  *1  '  S,1  +  b2°3*)3  n3Sal  +  blV2*l2  8"l) 


n„  ~  1 


n„  - 1 


n.  -1 


+  (b1e]n2a22  8a?  +  b^aj3  n3  5a2  +  Sa2^]dS  -  0 


or,  with  (1), 


rU U.  »\J0  \  , 

i/s/hrr-**,  +  r^*-9ldSi0 


W  u. 

»a]  1  »#2 


The  volume  of  the  reinforcement  can  be  expressed  as: 


V  =  C  y*  a  d  S,  where  C  is  a  positive  constant  and  a  =  a^  +  a2 
8V  =  C  y* 8 a  d5 

=  C /  (8a]  +  8a2)dS  (7) 

The  design  B  will  not  have  a  greater  reinforcement  weight  than  any  neighboring  design 
B  of  the  same  stiffness,  if  the  reinforcement  volume  of  B  is  not  greater  than  the  rein¬ 
forcement  volume  of  B .  Hence,  a  condition  has  to  be  found  for  which, 


or,  from  (7), 


2? 


5V  >  0 


y’(5a1  +  5a2)dS  >  0 


Then,  from  (6)  if  follows  fhaf  fhis  is  frue  if, 


dll.  dU, 


da,  da. 


-  C, 


(8) 


where  Cj  is  a  positive  constant,  the  optimality  constant. 


From  (8)  it  follows  that: 


Uj  =  C1  aj  +  constant 


where  the  constant  is  zero  since  Uj  =  0  for  =  0,  as  follows  from  (1).  Similarly, 
U2  =  C1  «2  -  Furthermore, 


2U  =  U1  +  U2  =  U,  (1  4  -£) 


U.  — 
1  al 


Thus,  the  optimality  condition  becomes: 


2U 

a 


u,  U2 

—  =  —  =  constant 
al  a2 


(9) 


Equation  (9)  is  a  sufficient  optimality  condition  in  the  sense  that  it  leads  to  a  reinforce¬ 
ment  that  is  not  heavier  than  any  other  reinforcement  of  the  same  stiffness. 


To  prove  the  necessity  of  the  optimality  condition,  the  external  work  Wg  is 
maximized  for  given  V.  This  leads  to  the  variational  expression: 
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/*  p  n.  n<>  no  i  3,i  1  \ 

8 J  ♦  b2.3.,3  ♦  b]V,2)  ♦  (b,.,.2  ♦  b2.3.2  ♦  b,.2.2  ) 

-  P.5.  -  AC  (a,  +  a2)]dS  =  0 
where  A  is  a  Lagrange  multiplier. 


The  Euler  equations  reduce  to: 


n„  - 1 


n2-l 


,v;!"  ♦  vav?" *  V.V?  (,ob) 


n„-  1 


n, "  1 


Combination  of  (10a)  and  (10b)  yields  (8). 


The  optimum  is  weak  and  relative.  Weak  because  of  the  assumption  made  for 
(5),  namely  that  5a.  be  very  small  [213.  Relative  for  two  reasons:  firstly,  because  of 
the  statical  indeterminacy,  as  is  discussed  in  Paragraph  7.5.  and  secondly,  because  only 
the  optimality  of  the  reinforcement  is  considered.  That  means  optimization  with  respect 
to  the  total  weight  of  fibers  and  matrix  together. 


8.3.  Example 

The  simplest  example  illustrating  the  potential  of  the  method  is  a  wedge  in  ten¬ 
sion  (see  Figure  1).  The  isotropic  case  has  been  discussed  in  [37]  and  the  anisotropic  case 
in  [38] .  Details  of  the  analysis  appear  in  the  Appendix.  Based  on  data  for  boron-epoxy, 
some  numerical  values  were  computed,  which  are  represented  graphically  (see  Figure  2, 
where  the  optimal  design  is  designated  by  A  and  the  unidirectional  one  by  B). 
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In  Figure  2,  Weight  it  the  weight  of  that  pert  of  the  wedge  thet  is  bounded  by 
±  .  .„d  r,  end  The  reciprocel  of  the  ..rein  energy  in  the,  per.  of  the  wedge  i.  . 

of  the  stiffness.  The  specific  stiffs  is  defined  ..  the  reciprocel  of  the  product 

of  the  strain  energy  times  the  weight. 

For  increasing  wedge  angle,  the  graph  for  the  stiffness  -hows  the  increasing 
superiority  of  the  optimal  design  over  the  unidirectional  design. 

In  the  graph  for  the  ratio  of  the  weight  of  the  optimal  design  to  the  weight  of  the 
unidirectional  daslgn,  a  distinct  decrease  become,  evident  at  approximate  10  degrees. 
This  can  be  explained  as  baing  due  to  the  fact  that  past  that  angle  the  cosine  starts  to 
difhr  substantially  from  1 ,  and  decease,  mom  rapidly.  Hr.  fiber  spacing,  bain,  a  func¬ 
tion  of  the  -cant,  therefor,  increase,  owe  rapidly  which  in  turn  cause,  the  weight  to 
incra...  at  a  decreasing  rate.  And  so  the  ratio  of  the  weight,  of  the  optimal  an  the 

unidirectional  design  decreases. 

Th.  fact  that  th.  optimal  deign  is  significantly  stiffer  a.  well  a.  lighter  than 
unidirectional  deign  is  illustrated  by  Hr.  curve  for  the  ratio  of  the  specific  stiff— 

9.  Fiber  LAying  Gadget  (FLAG} 

In  order  to  check  theoretical  reult.  with  experimental  data,  tet  specimen.  hav. 
ba  mad.  in  which  th.  fiber,  run  along  Pai.ch.rie  predetermined  by  the  theory.  A 
device,  called  FLAG,  was  mad.  to  lay  a  single  continuous  wire,  thread  or  fiber  along 
any  predetermined  curve.  In  this  section  this  FLAG  is  described.  Th.  FLAG  Itself  is 
mounted  on  the  pen  carriage  of  an  X-Y  recedes,  replacing  th.  pen  apembly.  Th.  X-Y 
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recorder  is  controlled  by  a  P DP-1 2  digital  computer  through  a  D/A  (digital  to  analog) 
converter. 

With  appropriate  modification!  the  FLAG  alio  can  be  used  for  curved  surface! 
of  almost  every  type.  Its  range  of  application,  therefore,  is  less  restricted  than  is  the 
case  for  fiber  winding  techniques. 

9.1 .  Description  of  the  FLAG 

Plumber's  delight,  pictured  in  Figures  3,  4,  5  and  6,  and  made  of  leftovers  and 
some  scrap  pieces  of  copper  tubes  and  brass,  is  FLAG,  an  experimental  Fiber  LAying 

Gadget. 

Details  of  its  construction  are  shown  in  Figure  5.  It  consists  of  a  wheel  (1) 
pressing  the  fiber  onto  a  sticky  layer  and  simultaneously  pulling  the  f  er  (3)  from  the 
reel  (2).  The  wheel  and  reel  are  rigidly  interconnected  by  the  tube  (4)  through  which 
the  fiber  is  guided  from  the  reel  to  the  wheai.  This  prevents  the  fiber  from  being 
twisted  when  the  wheel  follows  a  curved  traiectory  and  rotates  around  a  vertical  axis. 
The  wheel-reel  assembly  is  supported  by  two  aligned  bearings  (5),  allowing  for  unre¬ 
stricted  rotation  around  the  vertical  axis  of  the  aligned  bearings.  The  bearing  assembly 
is  fixed  to  one  end  of  a  curved  tube  (6)  by  means  of  two  joints  (7)  and  (8  which  have 
mutually  perpendicular  horizontal  axss.  Both  joints  are  rigid,  but  they  are  adjustable 
for  alignment  purposes.  The  other  end  of  the  curved  tube  is  fixed  to  a  horizontal  bit 
that  fits  in  the  pen  carriage  of  the  X-Y  recorder  (see  Figure  4).  Like  the  pen  assembly, 
the  bit  can  rotate  freely  around  its  horizontal  axis.  The  wheel  is  covered  with  toflnn 
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tape  to  prevent  it  from  sticking  to  Hie  layer.  To  avoid  practical  complications,  no  pro¬ 
vision  was  made  for  a  fiber-cutting  device. 

9.2.  Operation  of  the  FLAG 

The  FLAG,  mounted  on  the  pen  carriage  of  the  X-Y  recorder,  is  controlled  by 
a  program  in  the  PDP-12.  Coordinates  for  the  FLAG  are  calculated  in  digital  form, 
and  by  means  of  a  D/A  converter  these  digital  data  are  converted  into  voltages  which 
activate  the  X-Y  recorder.  The  potentiometers  of  the  recorder  provide  fine  adjustment, 
necessary  to  scale  the  computer  output  to  Hie  X  and  Y  ranges  of  travel  of  the  carriage. 

The  trajectory  of  the  wheel  is  the  fiber  trajectory  determined  through  theory. 

The  trajectory  of  the  FLAG  is  calculated  as  follows:  Let  b  be  the  dist  .nee  between 
the  vertical  axis  of  the  two  bearings  and  the  point  of  contact  between  wheel  and  layer 
(see  Figure  6);  Let  (x^ ,  y^ )  be  a  point  on  the  fiber  trajectory  and  (x^,  y^  the  corre¬ 
sponding  point  on  the  trajectory  of  the  vertical  axis.  Then, 

Xj  =  Xj  +  b  cos  TJ 

y2  =  Y\  +  b  cos1\ 
dy1 

where  t)  is  defined  by  tg  t)  =  -t — . 

*1 

Since  the  fiber  trajectory  is  known,  tj  can  be  calculated  for  any  point  (x^  ,  y  .. )  and 
hence  (Xj,  y^)  can  ba  determined.  This  computation  is  executed  in  the  PDP-12. 
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9,3.  Initial  Results 

In  order  to  determine  experimentally  the  effect  of  the  interspace  between  the 
fibers  on  the  elastic  constants  of  the  composite,  several  unidirectional  test  specimens 
have  to  be  made,  with  different  fiber  soacing.  To  this  end  a  program  for  the  'DP-12 
was  written  to  make  the  FLAG  lay  fibers  along  straight  parallel  lines.  The  distance 
between  the  lines,  i.e.  the  interspace  between  the  fibers,  can  be  adjusted  by  the  corre¬ 
sponding  potentiometer  on  the  X-Y  recorder. 

First  trial  runs  with  thread,  nylon  and  steel  wire  indicated  that  the  fiber  needed 
to  be  guided  almost  to  the  point  of  contact  between  whe'l  and  layer.  This  was  accom¬ 
plished  by  flattening  the  lower  end  of  the  fiber  guiding  tube  and  making  it  follow  the 
contour  of  the  wheel  (see  Figure  7).  Two  different  wheels  of  different  diameter  were 
used--the  bigger  one  for  the  stiff  wire  and  the  small  one  for  the  flexible  thread  and 
nylon. 

In  the  first  prototype  of  the  FLAG,  the  distance  between  the  vertical  axis  and 
the  point  of  contact  between  wheel  and  layer  was  about  twice  as  long  as  it  is  in  the 
current  version.  This  was  done  to  avoid  wheel  shimmy.  Consequently,  however,  the 
radii  of  tho  return  loops  at  the  end  of  each  trajectory  had  to  be  quite  large.  Therefore, 
the  distance  was  reduced  rather  arbitrarily  to  its  present  value.  Still  no  wheel  shimmy 
occurs. 


For  the  development  of  the  FLAG,  nylon  fiber  has  been  used  almost  exclusively 
because  of  easy  handling  combined  with  moderate  stiffness.  At  first,  masking  tape 
served  as  the  sticky  layer,  but  the  fiber  usually  sank  into  the  adhesive  layer  of  the  tape. 
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too  deeply  to  be  covered  sufficiently  by  the  resin  applied  afterwards.  Evaluation  of  other 
materials  finally  led  to  Booklon,  a  product  of  Boeketo  Plastics  N.  V.  in  the  Netherlands. 
Normally  used  to  cover  books,  maps,  etc.,  Booklon  is  a  transparent,  self-adhesive,  plas¬ 
tic  sheet  with  a  protective  covering  which  is  easily  peeled  off.  The  adhesive  layer  is  very 
thin  and  the  fibers  therefore  are  exposed  sufficiently  to  allow  satisfactory  wetting  during 
subsequent  application  of  resin.  The  protective  covering  performs  well  when  used  as  a 
releasing  top  layer.  Once  the  resin  is  cuted,  the  Booklon  base  layer  and  the  protective 
covering  can  be  peeled  off  easily.  The  resin  content  is  well  controlled  when  some  pres¬ 
sure  is  applied  during  curing,  since  the  fibers  act  as  spacers.  For  void  control,  the  method# 
commonly  employed  in  hand  lay-up  techniques  are  used. 

10.  Conclusion 

In  the  foregoing  a  method  has  been  discussed  to  obtain  relatively  optimal  rein¬ 
forcements  in  FRC  membranes.  It  is  based  on  the  observations  that: 

1.  In  FR^'s  the  direction  of  maximum  st'ffness  coincides  with  the  fiber 
direction,  and 

2.  Principal  trajectories  are  independent  of  fiber  reinforcement  as  long 
as  the  fibers  coincide  v^ith  the  trajectories. 

An  optimality  condition  is  derived  which  states  that  the  specific  strain  energy 
divided  by  the  specific  number  of  fibers  is  a  constant.  The  method  is  illustrated  by  the 
example  of  a  single-layer  wedge  in  tensior.  The  results  obtained  for  the  example  indi¬ 
cate  the  potential  of  the  method  proposed:  for  a  20-degree  wedge,  the  improvement  of 
the  specific  stiffness  can  be  in  the  order  of  50  percent.  Therefore,  it  will  be  worthwhile 
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to  investigate  more  complicated  cases  and  to  formulate  criteria  concerning  the 
necessity  of  a  locally  orthogonal  reinforcement  pattern  or  the  sufficiency  of  a 
locally  unidirectional  pattern.  It  will  also  be  worthwhile  to  try  to  extend  the  optim¬ 
ality  criterion  by  including  the  matrix  in  the  optimization. 

To  take  full  advantage  of  FRC's,  coupling  effects  have  to  be  included  because 
of  their  possibly  beneficial  consequences  (e.g.  with  regard  to  stability). 

These  points  illustrate  Hie  investigation  of  optimality  of  FRC  plates  is 
pertinent  to  both  the  theory  and  the  application  of  FRC's. 
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Figure  2.  Ratios  of  Specific  Stiffnesses, 
Stiffnesses  and  Weights. 
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(see  paragraph  9.3.) 


Appendix;  Analysis  of  the  Wedge 
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Solutions  for  an  isotropic  wedg«  and  an  anisotropic  wedgo,  respectively,  can  be 
found  in  [37],  p.  97  and  in  [38],  p.  87,  respectively. 

The  composite  is  assumed  to  be  an  epoxy  matrix  reinforced  with  boron  fiben. 
The  following  data  are  used: 

Ef  =  4.5  x  I06  kg/cm2 

E  =  0.035  x  I06  kg/cm2 

m 

Pj  =  2.35  g/cm^ 

,3 

P  =  1.15  g/cm 

m 

As  for  the  geometry  of  the  wedge  (see  figure  1)  /  r|  =3  and  0  <  a  <  15° 

The  optimal  solution  is  compared  with  the  solution  for  a  w@«%e  made  of 
unidirectional  tape,  oriented  in  the  x  -  direction,  with  the  following  constants; 

fiber  fraction:  65%  by  volume 

E^  =  2.8  x  10^  kg/cm2 

Ey  =  0.28  x  I06  kg/cm2 

V  =  0.25 

xy 

C  =  0.11  x  I06  kg/cm2 
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The  stress  distribution  In  an  isotropic  wedps  is  given  by  137]  : 

P  cos  9 

a  =  .  . . 

r  ( a  +  1/2  sin  2a  ) 


where  a  is  the  stress  in  the  r- direct! on  (figure  1):  all  other  stress  -  components 
are  zero.  The  wedge  is  assumed  to  be  of  unit  thickness. 

The  principal  trajectories  are  straight  lines  in  the  r  -  direction  and  concentric 
circles  in  the  9  -  direction. 

Fibers  are  laid  in  the  r  -  direction  only. 

From  the  optimality  condition,  paragraph  8.2.  (g),: 


U 


it  follows - -  C 

E  a 


(Al) 


where  E  =  E  -  modulus  in  the  r  -  direction 


Cj  -  constant 


C  = 


2  (a  +  1/2  sin  2a)2 


Tb.  number  of  flbm  Is 


(A  2) 


a3 


where  a  (  =  a  at 

f  (0)  =  function  of  0,  to  be  determined. 

For  E  the  rule  of  mixture*  holds,  hence  E  is  a  linear  function  of  the  fiber 
volume  fraction  and  consequently  also  a  linear  function  of  a,  the  number  of 
fibers,  as  is  shown  later  in  this  appendix. 


E  =  V  4  C4 


(A3) 


c3  -vn  ''ey 


=  fiber  volume  percentage  for  a  =  constant 


C.  =  E  ,  the  E  -  modulus  of  the  matrix 
*  m 

Ef  -  E  -  modulus  of  the  fiber 

Henceforth  will  be  omitted  since  E^<c  Ey  for  the  matefial  considered. 
Then,  from  (A  1 X  (A  2)  and  (A  3)  it  follows: 


C3a 


2  co«2  f 

cp~72 - 


(A  4) 


or  f  (0) 
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Substitution  of  (A  4)  into  (A  2)  gives* 


I  «G  •  Q  — 


cos  6 
r 


After  substitution  of  (A  5)  into  (A  3),  the  E  -  modulus  becomes 


assuming  E  «  E,. 

m  r 


(A  5) 


(A  6) 


The  constant  Cj  can  be  determined  from  the  boundary  condition  at  r  =  r^, 

=  0.  At  this  point  the  stress  is  maximum  and  therefore  the  number  of 
fibers  ought  to  be  maximum.  For  the  sake  of  simplicity  a]  is  assumed  to  be  1. 
(If  *|  >  1,  one  or  more  than  one  additional  layer  is  needed.) 


Expression  (A  5)  then  gives 


(A  7) 


\ 
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The  strain  follows  from  i  -  ;  thus: 


that  It,  tha  strain  In  tha  r  -direction  Is  constant 
Substituta  (A  7)  Into  (A  6),  raspactivaly,  to  obtain 


ri 

a  =  1  —  cos  0 

E  =  C3  ~  cos9 

(Still  under  tha  assumption  that  ij  =  1  at  6  =0  ). 

Tha  strain  energy  accumulated  In  tha  wedge  between  r  -  r}  and  r  =  r?  Is: 

UA  -  ffuiS- 

/»  r*n 

=  C7  /  a  rdr  de 


(A  8) 

(A  9) 
(A  10) 


Substitution  of  (A  7)  and  (A  9)  yields: 


Tha  relationship  between  fiber  fraction  and  number  nf  fibers  Is  found  as  follows: 


L*r  d  b*  the  fiber  diameter  and  Xf  the  distance  between  the  canters  of  two  neighbor¬ 
ing  fibers.  Then  a  *  d  /  since  a  Is  defined  as  the  number  c*  fibers  per 
unit  area  of  the  membrane  (the  unit  area  can  conveniently  be  chosen  as  being 
d  x  d). 

The  fiber  fraction  (volume  percentage)  per  unit  area  is: 

vf- . 

or  Vf  “  0.785  a 

The  total  weight  of  the  fibers  is 

Wf  =  //pf  Vf  rdr  d  0  = 

s  1.57  pf  r}  (r2-  r})  sina  (A  12) 

where  pf  *  specific  gravity  of  the  fiber  material 
The  total  weight  of  the  matrix  is 

wm  *  Sf»m  o  -  V  *  - 

=  Pm  ^r2*rl^  ^  (r2  +  rl  )Q  “  1.57  r^  slna]  (A  13) 

where  p^  =  specific  gravity  of  the  matrix  material 

(A  13)  added  to  (A  12)  yields  the  total  weight  (per  unit  thickness  d)  of  the  wedge 
between  r «  r^  and  r  = 


a7 


WA  =  W,  +  W  = 

A  f  m 

=  (r2  "  rj)  Cl^7  (  Pf-  p^)  rj  stna  +  Pm  +  r^a]  (A  14) 

For  a  wedge  of  anisotropic  material  a  solution  has  been  given  in  [38]. 

The  stress  distribution  is  given  by: 

_ 1  (A  cos  9  +  B  sin  0) 

"  -  7"  1 - LW— ■  -J ■  (A  15) 


where  a  ■  stress  in  the  r  -  direction;  all  other  stress  components  are  zero 
1  /  L  (0)  =  E  -  modulus  in  the  r  -  direction 

A,  B  =  constants  to  be  determined  from  the  equilibrium  conditions  in 
x  -  and  y-  direction: 


=  P 


r  d  0  =  0 


where  o*  rd0  -  a  cos  0  rd9 
a  rdO  =  a  sin  0  rd0 

y 


(A  16)  and  (A  17)  are  expanded: 


A 


coP  9 

TW 


do  + 


0  cos  0 


nur 


(A  16) 
(A  17) 


(A  18) 


d  0  =  P 


a8 


i 


+  a 


sin  9  cos  9 

TW— 


/♦•  2 

■T75y  d®  -  o 


L  (9)  ean  be  computer1  from  the  principal  constants  and  the  usual  transformation: 

L' 

y 


L(»)  -  ,I„29c«28  + 


By  numerical  Integration  for  a  set  of  a  values  ranging  from  1°  up  till  and  Including 
15°  the  corresponding  values  for  A  were  computed.  Because  of  the  symmetry  of 
the  geometry  with  respects  to  9  =  0,  and  since  L  (9)  Is  an  even  function  of  9 
while  sin  9  Is  not,  the  second  integral  In  (A  18)  and  the  first  one  in  (A  19)  are 
zero,  and  hence  B  =  0. 


Now,  an  expression  for  the  strain  energy  Is  derived: 


UB  ~ffBt  '*  <*' 
fr7  ra 

~l  /*E  tAco>9  +  Bsin®32  d9  = 

r,  j-a  a 

~  j  |  A  IaJ  E  cos2  9d9  +  B  j  E  sin  9  cos  9d9]  + 
f a  /*® 

+  B  CA  /  E  sin  9  cos  9d9  +  B  /  E  sin2  9  d9]  |  y 

Substitution  of  (AI8)  and  (AI9)  yields: 

UB  -  AP 

The  total  weight  of  the  anisotropic  wedge  (per  unit  thickness  d)  between  r  =  r. 


(A  19) 


(A  20) 


of 


and  r  =  r^  i« 

WB  "  5*3  (  M  +  Vm  "m*  (r22'rl2) 

( a  In  degrees) 

Substitution  of  the  data  for  the  geometry  and  for  boron-epoxy  unidirectional 
tape  as  given  at  the  beginning  of  this  appendix,  yields: 

WB  =  3^3  (0.65  X  2.35  +  0.35  x  1.15)  8  r}  2  x  I0“3  = 

=  0.27  x  10  3  x  rj  2  x  a  (kg) 

Since  in  the  point  (r  =  r^,  6  =  0)  the  stress  is  maximum  in  both  wedges,  with  the 
maximum  having  the  same  value  in  both  cases,  the  results  derived  for  the 
optimal  solution  have  to  be  adjusted.  The  value  of  Vj.  for  a  =  1  has  to  be  0.65 
instead  of  0.785.  Instead  of  (A  14),  the  following  expression  is  used  for  the 
calculations: 

WA  =  Vrl)  n'3  (W  Vlna  +  *m  Vi)el 

Substitution  of  the  data  for  the  geometry  and  for  boron  and  epoxy  gives: 

WA  =  r12  2C1.3  (2.35-  1.15)  sina  +  1.15  x  4  x  -^  ]  x  10'3  = 

=  [3.12  sina  +  0.16a]  x  r^  2  x  10  3  (kg) 

(  a  in  degrees) 


(A  21) 


(A  22) 


(A  23) 


(A  24) 


The  graphs  in  figure  2  were  calculated  using  (A  II),  (A  20),  (A  22)  and  (A  24) 
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